Long-Range Lepton Flavor Interactions and Neutrino Oscillations 
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Recent results from the MINOS accelerator neutrino experiment suggest a possible difference be- 
tween and Vp, disappearance oscillation parameters, which one may ascribe to a new long-distance 
potential acting on neutrinos. As a specific example, we consider a model with gauged B — L e — 2L T 
number that contains an extremely light new vector boson m z > < 10~ 18 eV and extraordinarily 
weak coupling a' < 10 -52 (or larger m Z ' if cosmology bounds on neutrino decay apply). In that 
case, differences between — > v T and v^ L — > v T oscillations can result from a long-range potential 
due to neutrons in the Earth and the Sun that distinguishes and v T on Earth, with a potential 
difference of ~ 6 x 1CP 14 eV, and changes sign for anti-neutrinos. We show that existing solar, 
reactor, accelerator, and atmospheric neutrino oscillation constraints can be largely accommodated 
for values of parameters that help explain the possible MINOS anomaly by this new physics, al- 
though there is some tension with atmospheric constraints. A long-range interaction, consistent 
with current bounds, could have very pronounced effects on atmospheric neutrino disappearance in 
the 15 - 40 GeV range that will be studied with the IceCube DeepCore array, currently in operation, 
and can have a significant effect on future high-precision long-baseline oscillation experiments that 
aim for ±1% sensitivity, in and disappearance, separately. Together, these experiments can 
extend the reach for new long-distance effects well beyond current bounds and test their relevance 
to the aforementioned MINOS anomaly. We also point out that long-range potentials originating 
from the Sun could lead to annual modulations of neutrino data at the percent level, due to the 
variation of the Earth-Sun distance. A similar phenomenology is shown to apply to other potential 
new gauge symmetries such as L — 3L T and B — 3L T . 



I. INTRODUCTION 

Neutrino flavor oscillation experiments have provided 
some of the most direct and robust indications of physics 
beyond the Standard Model (SM). Solar, atmospheric, 
reactor, and accelerator data all point to the conclu- 
sion that at least 2 active neutrinos have tiny but non- 
zero masses of up to order 0.1 eV, whose generation re- 
quires extending the SM. We refer the interested reader 
to Refs. U | for a review of the extensive literature 
on neutrino oscillation physics. Given the smallness of 
neutrino mass differences, even minute perturbations to 
the time evolution of flavor eigenstates, caused by fee- 
ble differences of interactions of neutrinos with back- 
ground sources, can produce measurable departures from 
vacuum oscillations. For example, these effects can be 
caused by the short-distance electroweak interactions of 
neutrinos with solar or terrestrial electrons, referred to 
as the Mikheev-Smirnov- Wolfenstein (MSW) effect §, @. 
The sensitivity of neutrino oscillations to such small ef- 
fects makes them a good probe of new physics that vio- 
lates u e -v^-v T universality ||. Hence it is interesting to 
look for unexpected effects in neutrino data. 

Recently, measurements at the MINOS experiment Q 
have resulted in different inferred values for differences of 
squared masses and mixing angles 



[where sin 2 (2^23) = 1-00 is the best fit value, while 
sin 2 (26» 23 ) > 0.91 at 90% confidence level] and 
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0.86 ±0.11 
(2) 

in and disappearance, respectively. The above MI- 
NOS results have revived some interest in long-range in- 
teractions (LRIs) that can cause disparities between 
neutrinos and anti-neutrinos. For other related works on 
the MINOS anomaly, see, for instance, Refs. jlO|] . 

The possibility of new long-range forces was discussed 
in the pioneering work of Rcf . [ p"l[ , and subsequently con- 
sidered as an alternative way to explain apparent CP vi- 
olating effects in K meson decays |l^, [l3|. Note that 
the disparity in the oscillation parameters for neutrinos 
and anti-neutrinos, as suggested by the MINOS results 
(|]) and (|j), can be ascribed to an apparent violation 
of CPT W%- However, in what follows we will assume 
that CPT is conserved in vacuo and consider the possi- 
bility that the MINOS result could be a hint of a new 
LRI. Eotvos-type ]l5| tests of gravity place stringent 
bounds on these interactions [pT] , constraining their "fine 
structure constant" a 1 < 10"^ (electron coupling) and 
a' < 10~ 47 (nucleon coupling) [I6| |r7|. This suggests 
an astrophysical source with a large number of parti- 
cles is needed, for sizable long-range effects. Long-range 
interactions in neutrino oscillations were considered in 
Refs. Jl8] ^l] ; see also Ref . |2^] . The long-range vector in- 
teraction yields equal and opposite potentials for leptons 
and anti-leptons. This can then manifest itself as a differ- 
ence in the properties of neutrinos and anti-neutrinos in 
terrestrial oscillation experiments, caused by the collec- 
tive effect of particles in the Sun and the Earth charged 
under a new U(l)' gauge symmetry. The correspond- 



ing effective fine structure constant must be extremely 
small, < 0(1O~ 49 — 10~ 47 ), as required by precision tests 
of gravity |23 . For comparison, note that the effective 
gravitational coupling between two protons is of order 
ot g ~ GNiTip ~ 10~ 38 , where Gn is Newton's constant 
and nip is the proton mass. Here, it is assumed that the 
associated Z' vector boson has a mass mz' < 1/Rav, 
where i? A u = 1 AU ~ 1.50 x 10 8 km - 10 18 cV" 1 is the 
mean Earth-Sun distance. (Later, we will limit our dis- 
cussion to values of m z > that are not far below 10 -18 eV, 
in order to exclude contributions from the rest of the 
Galaxy.) 

Before going further, we would like to make a few com- 
ments regarding the results (|]) and (||). First, the sug- 
gested MINOS anomaly is not at statistically significant 
levels, being at most a 2-sigma effect. In addition, the 
available atmospheric data from MINOS yield the best 
fit 24 (2-state mixing) 
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which does not support the above accelerator results, 
and, while also statistically limited, very mildly prefers 
an opposite sign for the effect. Taken together, the above 
considerations do not provide a strong case for invoking 
new physics. Nevertheless, we find the MINOS acceler- 
ator data sufficiently intriguing to motivate an exami- 
nation of the prospects for probing long-range leptonic 
forces at current and future experiments, as detailed be- 
low. 

In what follows, we will discuss the possibility of at- 
tributing the aforementioned MINOS anomaly to a LRI 
potential, generated by the neutrons in the Earth and 
the Sun. We will show that the existing bounds from 
neutrino oscillation data do not exclude such an inter- 
pretation. We use our approximate fit as a benchmark 
for potentially interesting values of parameters and es- 
timate the reach of current and future experiments for 
the new LRI. We find that the IceCube DeepCore array 
p5| , which is currently in operation, can provide an excel- 
lent probe of the benchmark model parameters and reach 
well-beyond them. We point out that long-range poten- 
tials generated by solar particles will inevitably lead to 
annual modulation of neutrino oscillations at Earth, due 
to the variation of the Earth-Sun distance. The large 
event sample expected at DeepCore seems sufficient to 
uncover a possible effect at the 1% level, statistically. 
Observation of such modulations can provide a distinct 
clue as to the solar contribution to the LRI and set a 
lower bound on its range. We will also consider long- 
baseline experiments, such as those envisioned for the 
Deep Underground Science and Engineering Laboratory 
(DUSEL) [E6|, to discover or constrain various effects of 
the LRI. We find that the expected capabilities of these 
experiments would allow them to probe the difference 
between the oscillation parameters of neutrinos and anti- 
neutrinos, induced by LRIs, which is a key signal for this 
type of new physics. Although we concentrate on the 
case of a potential generated by neutrons, our results and 



analysis carry over to other interesting scenarios where, 
for example, electrons or the total baryon number of the 
Sun and Earth may be responsible for the LRI. 

We will next briefly present the basic setup and formal- 
ism used in our work. Section III will contain our analysis 
and results. Our concluding remarks will be presented in 



Section [V 



II. FORMALISM 

Let us consider the addition of a general anomaly-free 
U(l)' gauge quantum number (for vectorial representa- 
tions) to the SM ||, H 

Q = a Q (B-L)+ai{L e -L ll )+a2{L e -L T )+a i {L ll -L T ), 

(4) 

where B and L arc baryon and lepton numbers, respec- 
tively, while Li, £ = e, fi,T are lepton flavor numbers, and 
dj, i — 0, 1, 2, 3, are arbitrary constants. For our primary 
example, we will set a\ = a-i = and for definiteness 
take ao = 03 = 1. However, any values of ao and 03 will 
lead to the same neutrino oscillation phenomenology for 
a fixed coupling between B — L and — L T . In this 
combination of quantum numbers, 

Q = (B - L) + (L„ - L r ) = B - L e - 2L Tl (5) 

(B — L) is associated with the source of the new potential, 
while (Lfj, — L T ) provides a contribution to the relevant 
neutrino oscillation — v r . Our choice for Q in Eq. (Q), 
as we will later argue, is less constrained by experiments 
than the previously studied L e — L^, £ = fx,T, cases. It 
also follows that the LRI potential due to B — L that we 
consider is generated by the total neutron number, since 
the contributions of electrons and protons cancel. 

Our charge assignment provides a simple way of achiev- 
ing the effective coupling in Ref. 0, where the micro- 
scopic origin of the requisite interactions is a mixing be- 
tween a Z' associated with — L T number and the Z 
boson of the SM. In principle, one could also imagine a 
mixing between two Z' states associated with, say, B — L 
and — L T , where an appropriate choice of mixing pa- 
rameters will yield the effective scenario adopted here. 
Given that our main purpose in this work is to eluci- 
date the relevant phenomenology, without reference to a 
particular underlying theoretical context, our choice of 
the gauged quantum number captures all the relevant 
key features for our analysis, while avoiding unnecessary 
complications. Note that as long as one of the anomaly- 
free quantum numbers is carried by a main constituent 
of solar or terrestrial matter, with the other lepton fla- 
vor number differences, one can build models that result 
in qualitatively similar effects. Indeed, we later consider 
models with gauged L — 3L T and B — 3L T that exhibit es- 
sentially the same phenomenology as that of B — L e —2L T , 
but with even smaller gauge couplings. 

The range of the interaction corresponding to charge 
Q is determined by the mass m z > of the force carrier Z' . 



2 



Since we are interested in the effect of a large number of 
particles, we assume that mz' ;$ 10~ 18 eV so that the 
neutrons both in the Earth and the Sun can contribute. 
We will not consider mz< <C 10~ 18 eV so that our as- 
sumed LRI does not extend far beyond the solar system 
and the contribution of the rest of the galaxy can be ig- 
nored p9| . The resulting potential felt by neutrinos on 
the Earth is then given by 
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using the estimated solar neutron fraction Y® = 1/7 (i.e. 
N®/N® ~ 6, where N® is the number of protons in 
the Sun) and Y® = 1/2 for the neutron fraction in the 
Earth. In Eq. (§), N® = 1.78 x 10 51 and N® = 1.70 x 
10 56 are numbers of neutrons in the Earth and the Sun, 
respectively, R® = 6.4 x 10 3 km is the Earth's radius, 
and Res is the variable distance of the Earth from the 
Sun. We note that i?ES attains its maximum value i?g S ~ 
1.52 x 10 s km at the aphelion (~ July 4) and its minimum 
value i?g S ~ 1.47 x 10 8 km at the perihelion (~ January 
4). 

The ratio of the potential V® at the Earth's surface 
from its neutrons to V® from solar neutrons is given by 
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Thus, the contribution of the Earth-generated potential 
is sub-dominant but not negligible. Note that if electrons 
are the source of the long-range potential, one can show 
that the solar contribution is roughly 24 times larger than 
that generated by the Earth BjJ . 

As discussed in Refs. || |l8|, p. 91, the survival proba- 
bility in the 2 flavor v^ — v T oscillation [for sm 2 (26*i3) ~ 0, 
Anij 2 ~ 0] is given by 



P w = l-sin 2 (20 23 )sin 2 



l 23 J 



(8) 



where 



Afh z 23 = Am' 23 J ~ cos{29 23 )} 2 + sm z (26 23 ) (9) 



and 
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Here, the symbols that are tilde-free denote vacuum 
quantities, and 
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with W T — Q T V n the potential energy for v T \ Q T = —2 
is the charge of u T , in our model. One can obtain the 



survival probability from the above expressions by 
f — > — f , and there is a degeneracy if both Am| 3 and 
cos(26* 23 ) change sign. Note that if sin (26*23) = 1 the 
formalism yields the same results for v and v. It should 
also be noted that the Am| 3 and sin 2 (26* 23 ) are energy- 
dependent for a 1 7^ and deviations from the vacuum 
values tend to increase with energy. 

The sin 2 (26' 23 ) and Am 23 measure the depth and lo- 
cation (Ey w Ato 23 L/27t) of the first oscillation mini- 
mum in the survival probability P w versus energy E v . 
With X = |f - cos(26» 23 )|, in the limit of X ~ (reso- 
nance condition), we have sin 2 (26 i 23 ) ~ 1 and Am 23 ~ 
Am 23 sin(26 , 23)- As X increases, sin 2 (2(? 23 ) decreases and 
Am| 3 increases. For a negligibly small cos(26' 23 ) where 
X ~ |f I = I - 2Q T V n E v /Aml 3 \, sin 2 (2023) decreases and 
Am^ 3 increases with E v for both v and v, for a given 
Am?, 3 . For a sizable cos(26 i 2 3 ), X may increase/decrease 
with E u depending on the relative sign of the two terms 
in X as long as |f| < | cos(26 ) 2 3 )|. This means devia- 
tions from the standard oscillations arc different for v 
and v. If a new potential (or in general X) is too large, 
sin 2 (26 l 2 3 ) ~ and the oscillation would be quenched. 



III. PHENOMENOLOGY 

In this section, we will examine the implications of new 
LRIs for current and future experiments. As a guide for 
our following discussion, we first derive, for the £7(1)' of 
Eq.([s]) where neutrons are responsible for distinguishing 

and v T , an approximate bound on a' based on the MI- 
NOS v ^ disappearance data (which dominate the statis- 
tics ||]). At 3-sigma, we roughly get a' < 5 x 10~ 52 , 
assuming cos(26>2 3 ) = [corresponding to the best fit 
value sin 2 (26 i 2 3 ) = 1 without new physics]. However, in 
order to address the disparity between the parameters of 
neutrinos and antineutrinos suggested by the MINOS re- 
sults (0) and (§), we will consider allowing cos (26*23) 7^ 
within the LRI scenario. We will next perform an approx- 
imate fit of the aforementioned MINOS results, obtained 
at a baseline of L = 735 km, within our reference model. 
Given the low statistical significance of the anti-neutrino 
results (~ 100 events ||), the fit is dominated by the neu- 
trino data points (~ 2000 events Q). For convenience, 
we employ the simplified MINOS data used in Ref. |L0| |, 
subtracting the neutral current background from the os- 
cillated signals. For the entire neutrino and anti- neutrino 
data set (23 bins) we find the best fit vacuum parameters: 
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2.4 x 10" 3 eV 2 ; sm 2 {29 23 ) = 0. 



1.0 x 10" 52 . 



(12) 



(13) 



with x 2 = 20.4 for 20 degrees of freedom. The above fit 
represents the preferred values of parameters in the pres- 
ence of the LRI, although we find that the goodness of 
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Ey (GeV) E r (GeV) E v (GeV) 

(a) (b) (c) 



FIG. 1: Survival probability P M(J for and z/ M with without (black dotted curve) and with the LRI (red solid curve for 
Up and green dashed curve for P M ), for Am^ = 2.4 x 10 -3 eV 2 , sin 2 

(2023) = 0.9, and a' = 1.0 x 10~ 52 . Typical values for the 
baselines have been chosen: (a) L = 2 x 6400 km (DeepCore), (b) L = 1300 km (DUSEL), and (c) L = 2 x 1300 km. The 
neutrino and anti-neutrino survival probabilities are different from each other in the presence of the LRI, since sin 2 (2^23) 7^ 1. 




E v (GeV) E r (GeV) E v (GeV) 

(a) (b) (c) 

FIG. 2: Same as Fig. |l| except for a' = 0.5 x 10~ 52 . 



fit is basically the same as the standard oscillations with 
no new physics; this was also the case in Ref. [[7 , where 
a fit but with a larger effective coupling was obtained. 
However, the parameters in Eqs. fll2| ) and (|u|) capture 
the implications of a new physics effect on the v ^ and 
data. Next, we will examine the implications of existing 
bounds for our fit to the MINOS results. 

First, let us consider the constraints from the solar and 
KamLAND data. The bound obtained in Refs. [§0[ || 
by comparing KamLAND and solar neutrino data leads 
roughly, for our neutron (rather than electron) based 
potential, to the increased bound a! < 6 x 2.5 x 
10- 53 /cos(26> 23 ) = 1.5 x Kr 52 /cos(26> 23 ) at 3 sigma; the 
factor of 6 comes from N® /iV® ~ 1/6, with TV® the num- 
ber of electrons in the Sun. For our value of cos(26>23) — 
0.3, the bound becomes a' < 5 x 10 -52 which is about the 
same as the rough MINOS bound given above and easily 
satisfied by our new physics scenario. To see how the 
quantity cos(2#23) enters into the bound with our choice 
of gauged quantum number, note that the solar neutrino 
oscillations can be described by two flavors: v e -R- v x . 
where v x = cos #23^1 — sin 023 ''•n assuming #13 — s- [the 
present bound is sin 2 (26>i3) < 0.15, at 90% confidence 
level A third eigenstate v y = sin 6*23^ + cos #23^1- 
decouplcs in this limit (for more details see Ref. pl|). 



With our choice of Q = B — L e — 2L T , we get 

(v e \H Lm \v e ) CX (l/elfieK) = -1 (14) 

{v x \HuaWx) k (v x \QxWx) = -2sin 2 9 23 , (15) 

where 7?lri is the contribution of the LRI to the Hamil- 
tonian. Since subtracting a matrix proportional to the 
identity in the evolution equation does not alter the os- 
cillations, we have effectively 











COS 26*23 



(16) 



for v e — v x oscillations, which yields the aforementioned 
suppression by cos(26 i 2 3). 

We now turn to the atmospheric constraints, which 
turn out to be the tightest. The effects of new physics 
on neutrino oscillations are often discussed in terms of 
coefficients ew pl"[ which parametrize the strength of 
the additional contributions in units of V2G]?n e , the 
standard MSW matter effect, with Gf the Fermi con- 
stant and n e the electron number density of the rele- 
vant medium. The analyses in Ref. |3l[] yields, at the 
95% confidence level, the upper bound e TT < 0.2. For 
atmospheric neutrinos traveling through the Earth, as- 
suming an average density of roughly 6 g cm -3 , corre- 
sponding to n e w 1.4 x 10 10 eV 3 , that constraint implies 
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FIG. 3: The black solid curves are numbers of upward moving atmospheric neutrino events per year at IceCube DeepCore as 
a function of energy, from Ref. |23] . The upper black curve corresponds to no oscillation and the lower black curve accounts 
for standard oscillations with sirr(2#23) = 1-00 and Am^ = 2.4 x 10 -3 eV 2 . The red dashed curve is our reproduction of 
the standard oscillation, with the assum ptio n of isotropic atmospheric neutrino flux with 0.687T < <j> < n, which is essentially 
indistinguishable from the result of Ref. |25| . 



W T = e T rV2G F n e < 4.6 x 10" 14 cV. Our MINOS fit cor- 
responds to W T ~ 5.6 x 10~ 14 eV. The value of W T from 
the LRI is however approximately constant throughout 
the Earth, due to the dominance of the solar contribu- 
tion, whereas the density of the Earth, sampled by neu- 
trinos traveling along its diameter, varies from roughly 
12 g cm -3 , corresponding to W T ~9x 10 -14 cV, in the 
core (which extends to ~ 3400 km from the center of the 
Earth) to 5 g cm -3 , corresponding to W T ~ 4 x 10 -14 eV, 
in the mantle (whose thickness is ~ 2900 km). Thus, the 
mean potential sampled by neutrinos traveling close to 
the diameter will be somewhere in between these val- 
ues, roughly W T ~ 6 x 10 -14 eV, somewhat larger than 
the value corresponding to our MINOS fit. Therefore, 
while our fit appears to have some tension with the atmo- 
spheric constraint, a more detailed analysis is called for, 
in order to pinpoint a precise constraint on the contribu- 
tion from the LRI we have assumed. Henceforth, we will 
adopt a' = 1.0 x 10 -52 as a plausible benchmark value for 
further exploration by experiments, where the potential 
MINOS anomaly suggests that new physics may appear. 
For our numerical illustrations, we will use vacuum pa- 



rameters At 



2.4 x 10" 3 eW and sin 2 



_-ai = 0.9, 

which are motivated by our fit results in Eq. ( |l2|) and 
are also reasonable in the standard scenario without new 
physics. In this way, we can focus on the effect of the LRI 
potential on neutrino oscillations, given specific vacuum 
parameters that could in principle be determined with 
high precision from other measurements. 

We also mention that the same type of study outlined 
above can easily be applied to other J7(l)' gauge groups 
with different quantum numbers and LRIs. Two interest- 
ing examples, gauged L — 3L T and B — 3L T , turn out to 
have essentially the same effect on MINOS as well as at- 
mospheric and long baseline neutrinos as our B—L e —2L T 
cxample, but with smaller values for a' of 1.3 x 10 -53 and 
1.1 x 10 -53 , respectively Here, the smaller values of a' 



are compensated by the larger numbers of electrons and 
baryons. 

With the completion of the DeepCore array, the Ice- 
Cube experiment will be able to probe atmospheric neu- 
trino oscillations for energies of order 10 GeV and above, 
well beyond the typical reach of the SupcrKamiokandc 
detector. The expected large statistics, of order 10 5 
events per year |^5| , make DeepCore an interesting probe 
of LRI using atmospheric neutrino data, which we now 
consider. 

Fig. |l|(a) shows the survival probability P MM ver- 
sus energy E v without (black dotted curve) and with the 
LRI, where the red solid and green dashed curves corre- 
spond to and v^, respectively. We have assumed the 
Earth diameter as the baseline, L = 2 x 6400 km, as a 
typical value relevant to the DeepCore array. Am| 3 = 
2.4 x 10" 3 eV 2 , sin 2 (2<9 23 ) = 0.9, and a' = 1.0 x 10~ 52 
have been chosen for illustration purposes. The sign of 
A77I23 nas Deen chosen according to our fit in Eq. jl^), 
and for definiteness we will choose cos(2#23) > through- 
out our analysis. We see that the LRI distinguishes neu- 
trinos and anti- neutrinos for sin 2 (2^23) 7^ L The de- 
viation of the LRI from the standard scenario becomes 
significant for E v > 15 GeV for neutrinos, whereas the ef- 
fect is quite large for anti-neutrinos over the entire range 
of energies considered here. 

Fig. [j](b) shows a similar plot for L = 1300 km, typical 
of the baseline for DUSEL experiments, which happens 
to be approximately 1/10 of the Earth diameter, with 
the other parameters as in panel (a). Deviations from 
standard oscillations do not become as significant as in 
case (a) for the same value of E v /L, due to the additional 
energy dependence of the LRI parameters in Eqs. ([)]) and 
( |l0| ) . The difference between the and signal is larger 
than their individual deviation from standard oscillations 
and can be a potentially distinct signal. 

Fig. |l](c) shows the effect of increasing the baseline by 
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FIG. 4: Atmospheric neutrinos per year (per 3 GeV) with E v at the IceCube DeepCore, using the same assumptions applied 
for Fig. ^[ Shown are the unoscillated case (top black dashed curves) and the case of no new physics (bottom thin black solid 
curves), as well as the cases a 1 = 1.0, 0.5, 0.1 x 10~ 52 corresponding to thick solid, dashed, and dotted curves, respectively. The 
vacuum parameters are sm 2 (26<23) = 0.9 with (a) Aml 3 = 2.4 x 10~ 3 eV 2 and (b) Am 2 3 = -2.4 x 10" 3 eV 2 . 



a factor of two, L = 2 x 1300 km, compared to the case 
presented in Fig. Mb). For the given energy range, as 
expected from Eq. (g), more oscillations will take place. 
We also see that the effect of the LRI on oscillations is 
much more prominent, but at about twice the value of 
E v compared to that in Fig. [|(b). It should be noted 
that merely increasing the baseline by a factor of two 
compared to the case in Fig. [l](b) would result in a re- 
duction by 1/4 in the flux and hence the event rate, all 
other factors being equal. 

The results in Fig. [I] suggest that it might be easier 
to see the LRI effect in the DeepCore experiment than 
at a long-baseline experiment. On the other hand, long- 
baseline experiments, unlike the DeepCore, can in princi- 
ple detect the asymmetry in the and oscillations for 
sin 2 (26*23) 7^ 1) which is a key feature of the LRI assumed 
here. 

Fig. ^ shows the same qualitative features for a smaller 
coupling a' = 0.5 x 10~ 52 ; as expected, the effect is less 
pronounced. In any event, the plots suggest that for 
E v > 15 GeV DeepCore can be quite sensitive to the 
new physics. Note that the effect on and can be 
interchanged by changing the sign of Ato 2 ^. 

To estimate the reach for LRI effects in neutrino oscil- 
lations at the DeepCore array, we will only consider the 
upward moving flux of atmospheric neutrinos. We use 
the no-oscillation number of events per year as a func- 
tion of energy, provided by the IceCube collaboration in 
Ref. (where these events are referred to as "verti- 
cal"). This event rate is represented by the top black 
solid curve (adapted from Ref. p5[ ) in Fig. ||. The lower 
black solid curve in this figure is the expected number 
of upward moving events for standard oscillations with 
sin 2 (20 23 ) = 1-00 and Am^ 3 = 2.4 x 10~ 3 eV 2 , also pro- 
vided by the IceCube collaboration ^5|. In order to re- 
produce this curve, we assume an isotropic flux of atmo- 
spheric neutrinos at each point around the Earth, and 



vary the zenith angle interval <po < § < 7r over which 
we calculate the event rate for upward moving neutrinos 
with oscillations, for the same values of sin 2 (26*23) and 
Am 2 ^ as above. The zenith angle <f> = 7t (tt/2) corre- 
sponds to the neutrinos moving vertically upward (hor- 
izontally) towards the detector. We have also assumed 
the ratio of to events with no oscillations to be 
2:1, reflecting a good approximation for the ratio of the 
cross sections Q . The red dashed curve in Fig. || is our 
reproduction of the oscillation with O = O.6871", which 
agrees very well with the solid curve (from Ref. |25j) and 
is basically indistinguishable from it in the plot. Hence, 
we will use the above approximations, with 4>q = 0.687T, 
to estimate the reach for new LRI effects at the DeepCore 
array, in the following. 

Fig. ||(a) illustrates the effect of the LRI on the num- 
ber of events (for a one-year run and per 3 GeV en- 
ergy bins) in and disappearance experiments. The 
top (dashed) curve corresponds to no oscillations p5| . 
The thick (red) solid, dashed, and dotted curves are for 
oi = 1.0, 0.5, 0.1 x 10 -52 , respectively and a' = (no new 
physics) is represented by the thin (black) solid curve at 
the bottom. The vacuum parameters are sin 2 (26*23) = 0.9 
and Am 2 ^ = 2.4 x 10~ 3 eV 2 , for all cases, to illustrate 
the effect of the LRI on neutrino oscillations. We see 
that except for the smallest value of a 1 the other cases 
are very distinct from the no new physics case, for 15 
< E u (GeV) < 40. Over this range of energies, the 
no-new-physics case yields roughly 3680 events and for 
a' = 0.1 x 10~ 52 the number of events changes by about 
100. Hence, statistically, a few-percent-level measure- 
ment, which might require 3 years of data, could in prin- 
ciple reach one order of magnitude below our benchmark 
value of a' = 1.0 x 10~ 52 , at the 3 sigma level. 

Fig. |4|(b) is the same as Fig. |I](a) except that Am| 3 = 
—2.4 x 10 -3 eV 2 , interchanging v and i>, is chosen. The 
thick (blue) solid, dashed, and dotted lines correspond to 
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FIG. 5: Percentage of annual modulation \(N a — N p )/(N a + N p )\ in atmospheric neutrino oscillation at DeepCore versus the 
LRI coupling a', for 15 < E v (GeV) < 30 around aphelion and perihelion. For each season, 120 days have been included. Other 
parameter values are the same as those of Fig. 0(a). 



the same values of a 1 as in Fig. [|(a) and the thin black 
solid curve at the bottom represents no new physics. We 
see that the effect of the LRI is now distinct for all values 
of a'. For 15 < E v (GcV) < 40, the number of events for 
the new-physics case with a' — 0.1 x 10~ 52 differs from 
that of the no-new-physics case by about 220. We see 
that even a several-percent-level measurement could in 
principle reach one order of magnitude below our bench- 
mark value of a! = 1.0 x 10~ 52 . Given the large statisti- 
cal samples expected at DeepCore, our estimates suggest 
that values of a' around one order of magnitude below 
that of our benchmark a' could potentially be probed by 
this experiment. 

Fig. H shows our estimate for the size of the an- 
nual modulation of atmospheric neutrino oscillations at 
DeepCore, as a function of the LRI coupling a'. Here, 
the vertical axis is \(N a — N p )/(N a + N p )\, where N a 
and N p are the numbers of events associated with aphe- 
lion and perihelion, respectively. An energy cut of 
15 < E v (GcV) < 30 has been implemented and 120 
days have been included around each apsis. For our es- 
timates, we have simulated the variation in Res by a 
sinusoidal function. This approximation of Earth's true 
Keplerian orbit captures the main effect we would like to 
illustrate, at the level of our analysis. The same param- 
eter values as in Fig. |(a), i.e. Am| 3 = 2.4 x 10~ 3 eV 2 
and sin 2 (2^23) = 0.9 are assumed. The total number 
of events (N a + N p ) per year for a' = 1.0 x 10~ 52 and 
a' = 0.5 x 10~ 52 arc about 2700 and 1400, respectively, 
and the seasonal modulations are near 1.2% for the for- 
mer and 0.8% for the latter. Our estimates then sug- 
gest that, depending on the value of a', 3-10 years of 
data could yield the necessary statistics to measure such 
levels of modulation. We have not accounted for atmo- 
spheric neutrino flux uncertainties, which can be as much 
as 10 — 15%. However, the (v^ + Vp)j{u e + v e ) flux ratio, 
which is proportional to our result, can be known much 
more precisely and would have uncertainties at the 1 — 2% 
level p3|. The large amount of statistics expected at 



DeepCore makes per-cent level measurements a realistic 
possibility |$4) . We may also expect that a more detailed 
and optimized analysis of the real data, using the pre- 
dicted time evolution of a stable flux ratio, could allow 
for a larger observed effect. Hence, our estimate suggests 
that DeepCore could be sensitive to annual modulations, 
when the solar source particles dominate the LRI poten- 
tial, as we have assumed in this work. 

Fig. ^| shows the predictions for the number of events 
per 0.125 GeV energy bins, over a 5-year run, for (a) 
and (b) long-basclinc disappearance experiments, e.g. 
at DUSEL, with L = 1300 km. We have taken the un- 
oscillated beam profile (dashed black curves at the top) 
from Ref. ||q] , corresponding to a 200 kt water Cerenkov 
detector. The vacuum parameters are sin 2 (26*23) = 0.9 
and Amj 3 = 2.4 x 10~ 3 cV 2 . The thin black solid 
curves correspond to no new physics. The thick red 
solid, dashed, and dotted curves again correspond to 
a' = 1.0,0.5,0.1 x 10~ 52 . We see that the largest two 
values of a' yield predictions that are distinct from no 
new physics, for 2 < E v (GeV) < 3 in both neutrino 
and anti-neutrino cases. We see that the LRI leads to 
distinct effects for and v^, which is a key signature 
of this new physics. However, the smallest value of a' 
does not yield a visible effect on these plots and may be 
difficult to reach at these experiments. 

Fig. [?] contains the same information, except for 
Am| 3 = —2.4 x 10~ 3 eV 2 , with the same conventions 
(new physics contributions are now displayed with blue 
thick lines). The same qualitative features as in the 
previous case with Am 2 3 > are present and reaching 
a' ~ 0.1 x 10 -52 seems to be difficult here as well. 

Fig. U shows the seasonal modulation for a 5-year 
run of the DUSEL as a function of the LRI coupling 
a' for i>n (red solid curve) and (green dashed curve). 
We take 2 < E v (GcV) < 3, with 180 days around 
each apsis. The same parameter values as in Fig. [|, 
sin 2 (26 l 23) = 0.9, Am^ 3 = 2.4 x 10~ 3 cV 2 , arc assumed. 
For a 1 = l.Ox 10~ 52 , the total number of events (N a +N p ) 
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FIG. 6: The number of (a) neutrino and (b) anti-neutrino events for a 5-year run (per 0.125 GeV) versus E v , in a long-baseline 
experiment with L — 1300 km (DUSEL). The unoscillated case (top black dashed curves) and the case of no new physics (thin 
black solid curves) are displayed, as well as the cases with a' = 1.0,0.5,0.1 x 10~ j2 corresponding to thick solid, dashed, and 
dotted curves, respectively. The vacuum parameters are sin 2 (2^23) = 0.9 and Am^ = 2.4 x 10 -3 eV 2 . 
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FIG. 7: Same as Fig. g except for Am| 3 = -2.4 x 10~ 3 eV 2 



are about 1100 (v^) and 2100 (v^) with the modulation 
at 0.5% and 0.4%, respectively. For a' = 0.5 x 10" 52 , 
the modulation is at 0.2%, for the total number of events 
about 1500 (i/ M ) and 1700 respectively. We assume 
the beam flux is constant over time. While our analysis 
is only a rough estimate, we may conclude that observing 
the modulation at the per-cent level would require exper- 
iments (beams and detectors) with somewhat more en- 
hanced capabilities, compared to those assumed for this 
analysis. 



IV. SUMMARY AND CONCLUSIONS 

In this work, we examined the effect of long-range in- 
teractions on neutrino oscillation experiments, motivated 
by the recent accelerator data from MINOS. These data 
are not conclusive, but could be suggesting that the oscil- 
lation parameters for neutrinos and anti-neutrinos may 



be distinct. Such an effect, if confirmed with more data 
in the future, could in principle be caused by long range 
interactions coupled to neutrinos (?], m M. As an il- 
lustrative example, we considered a U(l)' model with an 
ultra light gauge boson coupled to (B — L) + (L^ — L T ) = 
B — L e — 2L T that captures the key aspects of the req- 
uisite phenomenology. Such an effective interaction can 
also arise in other ways, for example through gauge boson 
mixing from two separate sectors Q] . The main required 
features are an interaction with a range of order 1 AU 
and coupling to both stellar matter and neutrinos, with 
the resulting potential characterized by a fine structure 
constant a' < 10 -52 . Alternatively, for a shorter range 
interaction, mz' > 10 -16 eV, the Earth itself may be the 
dominant source of terrestrial neutrino potential differ- 
ences. We pointed out that when the Sun is the dominant 
source of the long-range potential, the effect on neutrino 
oscillations at the Earth will exhibit annual modulations, 
due to the variable distance between the Earth and the 
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FIG. 8: Percentage of annual modulation \(N a — N p )/ (N a + N p ) | for a 5-year run versus the LRI coupling a', for 2 < E v (GeV) < 
3 around aphelion and perihelion. In each season, 180 days have been included, for (red solid curve) and (green dashed 
curve). The same parameter values as Fig. [] are used. 



Sun. 

We performed an approximate fit to the MINOS data 
within our reference model that gave a qualitative de- 
scription of the data. Our fit results accommodate the 
current experimental bounds on non-standard contribu- 
tions to neutrino oscillation data, although they show 
some tension with atmospheric bounds. However, for 
benchmark parameters motivated by our MINOS fit, we 
show that ongoing and future experiments could detect 
large effects due to the LRI potential, or else significantly 
further constrain such new physics. In particular, we esti- 
mated that the currently operational IccCube DecpCorc 
array can reach well beyond, by about an order of mag- 
nitude, the benchmark parameter space suggested by our 
MINOS fit, with about one year of atmospheric neutrino 
data at a typical baseline given by the size of the Earth, 
for 15 < E v (GeV) < 40. In addition, the large statistics 
afforded by DeepCore seem sufficient to detect a per-cent 
level modulation of neutrino oscillations with 3-10 years 
of data, providing key evidence for the solar source of the 
long-range potential. 

While DeepCore is only sensitive to the sum of neutrino 
and anti-neutrino events, they can be separately probed 
at future precision long-baseline experiments, such as 
those at a future DUSEL facility. With typical as- 
sumptions about the capabilities of such experiments, 
we showed that the values of parameters motivated by 
the MINOS results will be well covered, with about 5 
years of data. These experiments can in principle un- 
cover an asymmetry in the properties of neutrinos and 
anti-neutrinos, which is a key feature of the long-range 
interactions we have considered. Our simple estimates 
suggest the annual modulation of the data may not be 
easily accessible in these experiments unless their capa- 
bilities are somewhat more enhanced compared to our 
reference values. 

For /7(l)'s with gauged L — 3L T and B — 3L T , we found 
essentially the same phenomenology but with couplings 
a' about an order of magnitude smaller because of the 



larger sources of electrons and protons in the Sun, com- 
pared to neutrons. In all cases, DeepCore and future long 
baseline neutrino oscillation experiments are expected to 
push the a' sensitivity about an order of magnitude be- 
yond our benchmark values. 

Before closing, we would like to comment on some re- 
lated possibilities that could be of interest in the context 
of the new physics considered here. First of all, in our dis- 
cussion of annual modulation, we mainly considered fu- 
ture data. However, the Super-Kamiokande atmospheric 
data could already contain the annual modulation sig- 
nal in the ratio of the muon and electron neutrinos and 
a dedicated analysis could be of interest in view of the 
MINOS anomaly. Also, while we largely focused on a 
particular type of long-range interaction, interesting mo- 
tivation can be found in other possibilities. For example, 
the B — xiLi type of model, where the lepton sector has 
flavor-dependent couplings, can contain residual discrete 
symmetries that stabilize the proton and a dark matter 
candidate |2^]. This could provide possible new connec- 
tions between neutrino experiments and other areas of 
particle physics. It is conceivable that dark matter is 
charged under a new long-range force |36|, via a gen- 

or as in Rcf. [B8|. 



cralization of the setup in Ref. 37 



Then, new effects could perhaps arise from dark mat- 
ter trapped in the Sun or the Earth, in such scenarios. 
(While the modulation observed at the DAM A/LIBRA 
experiment ]39| ] might be explained by the motion of the 
Earth through the dark matter halo of the Galaxy, we 
find the relative proximity of the dates of the apses to 
the extrema observed by this experiment intriguing. A 
long-range potential from solar sources that also acts on 
dark matter could in principle affect the interpretation 
and implications of these measurements.) 

Hence, we conclude that current and future atmo- 
spheric and long-baseline neutrino oscillation experi- 
ments will provide an opportunity to probe long-range 
interactions whose feebleness generally excludes other ex- 
perimental search avenues. Together, such experiments, 
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like those at the DeepCore array and DUSEL, can test 
the relevance of these interactions to the suggested MI- 
NOS anomaly, or place more stringent bounds on their 
parameters. Detection of a new long-range force would 
be an important discovery with an immense impact on 
our view of fundamental physics. 
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